Ultrafast Gates for Single Atomic Qubits 
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Abstract 

We demonstrate single qubit operations on a trapped atom hyperflne qubit using a single ultrafast 
pulse from a mode-locked laser. We shape the pulse from the laser and perform a n rotation of 
the qubit in less than 50 ps with a population transfer exceeding 99% and negligible effects from 
spontaneous emission or ac Stark shifts. The gate time is significantly shorter than the period 
of atomic motion in the trap (^R a bi/^trap > 10 4 ), demonstrating that this interaction takes place 
deep within the strong excitation regime. 



PACS numbers: 03.67.Lx, 37.10.Vz, 37.25.+k 



Quantum information processing requires the ability to perform operations in an amount 
of time shorter than the coherence time of the qubit. The ratio between coherence time 
and gate duration can be increased by improving coherence times or developing faster gates. 
One experimental extreme is found in atomic systems such as trapped ions, which can have 
coherence times of many minutes l|, |2) . However, the precise spectral resolution needed to 
address narrow transitions and available excitation power have typically limited gate times to 
around a microsecond. In many condensed matter systems, on the other hand, single qubit 
operations must be performed on picosecond timescales [3J45J owing to the short coherence 
time of the qubit (T 2 * < 10 ns in |3j-l5|). Here we realize ultrafast operations in an atomic 
system with very slow decoherence rates. 

Aside from increasing the clock speed to decoherence rate ratio, there are also specific 
applications that would benefit immediately from fast single qubit gates, such as noise- 
reduction techniques [6], entanglement of material qubits with photon time-bin qubits [7|, 
and probabilistic gates with repetition rates limited by qubit rotations [8|. In particular, the 
ability to apply fast spin-dependent momentum kicks to trapped ions in the strong excitation 
regime |9| is a critical ingredient for fast sideband cooling [10] and fast entanglement of 



multiple atomic ion qubits ll|, [l2| . Such entangling gates can be performed faster than a trap 



oscillation period, in contrast to motional gates using spectroscopically resolved sidebands. 

We previously reported the implementation of an optical frequency comb to perform qubit 
operations 13]. In the low excitation regime, each pulse alters the atomic state by a tiny 
amount and the effect of the pulse train is manifested through careful coherent accumulation 
of spectral density into narrow comb teeth, resulting in gate times of order 100 /xs. In this 
Letter, we realize the strong excitation regime (^Rabi 3> zArap) [9] and demonstrate fast 
control of an atomic qubit by driving stimulated Raman transitions with a single pulse from 
a picosecond mode-locked laser. In order to implement a fast 7r-pulse, we perform simple 
pulse shaping in the form of a beam splitter and delay line [14J to fully rotate the qubit in less 
than 50 ps with population transfer exceeding 99%. Since the atomic qubit is well isolated 
from its environment, the time required to perform such a fast gate is only a small fraction 
(< 10~ 8 ) of the measured coherence time of our qubit. By setting the delay to zero in a 
counter-propagating geometry, the transition probability becomes sensitive to the atomic 
motion, which is the first step toward implementing fast cooling [10| and fast entangling 



gates between multiple ions [H, 
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FIG. 1. Experimental schematic showing (a) the energy level diagram of 171 Yb + with relevant 
states (not to scale), (b) the envelope of the electric field autocorrelation from the mode-locked 
laser and (c) the pulsed laser beam path with optional (dashed line) pulse-shaping beam path. 

Our experimental apparatus is shown in Fig. QJc). A m Yb + ion is trapped in a linear 
radio frequency Paul trap. The irtp = "clock states" of the two 2 Si/2 hyperfine levels act 
as the basis states of our qubit separated in frequency by z/hf = 12.642815 GHz. Doppler 
cooling, state preparation, and state detection are accomplished on the 2 Pi/2 -H- 2 Si/ 2 "Dl 
line" at 370 nm as described in Ref. 15j. The mode-locked laser for Raman transitions is a 



frequency-tripled yttrium vanadate laser operating at a repetition rate of /, 



rep. 



121 MHz 



with an average power of 4 W at 355 nm [16[. Single pulses are extracted by a Pockels cell 
pulse picker [17J and focused onto the ion. 

Fig. Ufa) shows the relevant energy levels of m Yb + . The 2 Si/2 qubit states |0) and |1) 
are coupled through excited 2 P states via far-detuned light that is polarized to drive either 
o~ + or o~ transitions, as Raman transitions from 7r light are forbidden by selection rules. 
The a ± Raman transition (two-photon) Rabi frequency is [18l Il9( 
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where g is the resonant one-photon Rabi frequency of the 2 P3/2(F = 2, m F = 2) -H- 2 Si/ 2 (-F = 
1, rrip = 1) "D2 line" cycling transition, A D1 is the detuning of the light above the 2 Pi/2 



state, Aps — Adi is the detuning of the light below the 2 P3/2 state, and Aps ~ 100 THz is 
the fine structure splitting. The ac Stark shifts of the qubit states from a ± light are jig, ll9j 

60 = 12 { AdT " A FS -A D J (2) 

6 = i_ ( 1 2 \ 

1 12VA D1 + z/ HF A fs -A d1 -z/hf/' 
where we have neglected the hyperfine splitting of the 2 P3/2 state. 

Eqs. HHshow that if the Raman laser is tuned between the Dl and D2 lines (0 < Adi < 
A FS ) the stimulated Raman transition amplitudes due to the 2 Pi/2 and 2 P3/2 couplings add 
constructively while the ac Stark shift contributions interfere destructively. The Stark shifts 
of the qubit states each cross zero near an optimal wavelength of A D P X ' = Aps/3, which for 
Yb + corresponds to a Raman laser wavelength of A opt . = 355 nm. The differential ac Stark 
shift of the qubit states does not exactly cancel, but has a local minimum near A op t. , reaching 
a value of 5o — S\ ~ fi ,i x 3z/hf/2Afs, corresponding to 2 x 10~ 4 Oo,i for Yb + . We measure 
a differential Stark shift of 1.1(5) x 1CT 4 fi ,i through microwave Ramsey spectroscopy with 
linearly polarized 355 nm light. 

The spontaneous emission rate can be estimated by calculating the excited state popu- 
lations during the Raman transition to give jig, [l9j 

rspon = 7 i [ A|7 + (A FS -A D1 ) 2 ) ' (4) 



where 7 ~ 2n x 20 MHz is the spontaneous emission rate from the 2 P states 20j. The 
probability of a spontaneous emission event during a 7r-pulse near A op t. can be estimated as 
"Pspon.Tr ~ 37/2 Aps- The local minima in r spon and P S pon,7r lie close to A op t. for Yb + at 349 
nm and 352 nm, respectively. At 355 nm we estimate V spon ^ n < 10~ 5 . 

In order to operate in this low Stark shift, low spontaneous emission regime, sufficient 
optical power at A op t. must be provided to make f2 ,i appreciable. For Yb + , frequency tripled 
mode-locked Nd YAG and vanadate lasers are available that provide many Watts at a center 



wavelength of 355 nm 21]. 



The two-photon Rabi frequency in Eq. [T]is time- dependent due to the shape of the pulse. 
As long as the pulse bandwidth is small compared to the single photon detunings A, the 
adiabatic elimination of the excited states remains valid and we can treat the qubit as a 
two-level system with a time- dependent coupling: 

H c s/h = — o z — a x . (5) 



Numerical solutions to the Schrodinger equation with this Hamiltonian can be obtained 
for a general time- dependent Q(t). However, there is an analytic solution due to Rosen and 



Zener for a coupling with a hyperbolic secant time-dependence [22J, which is the electric field 
envelope expected from picosecond mode- locked pulses [23(. Fig. d^b) shows the envelope 
of an electric field autocorrelation of a single pulse as measured with a scanning Mach- 
Zender interferometer and a fast photodiode, which is consistent with the autocorrelation of 
a hyperbolic secant with T pu i se = 14.8 ps having a linear pulse chirp of 8 x 10~ 3 ps -2 . Eq. 
[TJ however, shows that the two-photon Rabi frequency is proportional to the square of the 
electric field envelope. Nonetheless, sech 2 and sech are sufficiently similar that numerical 
solutions to the Schrodinger equation with O(t) oc sech 2 [2£/T] match the analytic solutions 
with Q(t) oc sech[7rt/T] to within a few percent for all simulations shown here. The Rosen- 
Zener transition probability from an initial state |0) to |1) is [22| 

V ^i = sin 2 (7ra) sech 2 [7rz/H F T pu i se ] (6) 

for a two-photon Rabi frequency of Qi o(t) = T u sech L 71 "* 1 where u is the pulse area. Eq. 

-* pulse ■*■ pulse 

|6] indicates that in order to create a high-fidelity 7r-pulse, the pulse duration must be many 
times shorter than the hyperfine period. 

Fig. [2]^a) shows the measured transition probability as a function of pulse energy (pro- 
portional to u in Eq. [6] and monitored with a fast photodiode). Both plots in Fig. [2] were 
taken with a single beam path, as represented by the solid line in Fig. Q^c). The maximum 
transition probability of 72% corresponds (see Eq. [6]) to a pulse duration of T pulsc = 14.8 ps. 

We also studied the dynamics of a train of identical pulses (with pulse energy denoted by 
the arrow in Fig. |2]^a)) at the laser repetition rate, shown in Fig. |2]^b). The solid curve is 
an analytic solution including free evolution with a fitted Gaussian decay in contrast. The 
fitted pulse duration is T pu i se = 14 ps, in good agreement with the fitted pulse width from 
Fig. [27a). 

It is clear from Eq. [6] and Fig. |S(a) that for T pu i se = 14.8 ps, even with unlimited pulse 
energy the population transfer probability cannot exceed 72% without changing the pulse 
shape. In order to execute a full 7r-pulse, we introduced a sim ple p ulse-shaping element into 



the beam path in the form of a beamsplitter and a delay line JjJ, resulting in two counter- 
propagating pulses as shown in Fig. QJc). By setting the energy of each half of the pulse 
to transfer 50% of the population, we are able to perform a fast Ramsey experiment with a 
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FIG. 2. (a) Qubit transition probability from a single pulse vs. pulse energy. The solid curve 
is a fit to Eq. [6] and the peak corresponds to a pulse energy of 12 ± 2 nJ. (b) Qubit transition 
probability vs. the number of identical (~ 8 nJ) pulses at the laser repetition rate. The solid curve 
is a fit to the Rosen-Zener solution with a Gaussian contrast decay. The single pulse energy for 
(b) is shown with an arrow in (a). 



variable free evolution time between two tt/2 Raman pulses. A 7r-pulse then corresponds to 
the top of a Ramsey fringe. 

The results of the pulse shape delay scan are shown in Fig. 02(a)- (c). We show results 
from experiments with three different polarization orientations for the two 7r/2-pulses. In 
(a), both beams drive only a + transitions. The Ramsey fringes are first maximized for a 
delay of 72 ps, corresponding to a net 7r-pulse after 80 ps. In the frequency domain, the 
shaped pulse spectrum in this case is a crude frequency comb with sinusoidal teeth separated 
by 14 GHz. The reason this separation is not exactly 12.6 GHz = z/hf is the additional a z 
rotation introduced by the Rosen-Zener solution dynamics. In the limit of infinitesimally 
short pulses, the comb spacing would converge to ^hf- Fig. MJo) shows the same experiment 
with one beam driving a + transitions and the other driving a~. The two central peaks each 
represent a ~ 40 ps net 7r-pulse, which is shorter than the 80 ps pure a + 7r-pulse due to 
the fact that the Rabi frequency for er + and o~ transitions have opposite sign and therefore 
rotate the Bloch vector about different axes. Fig. [31(c) shows the same experiment with 
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FIG. 3. Ramsey fringes for three different polarization orientations for each half of the shaped 
pulse obtained by fitting ensemble histograms of photon counts. The polarization configurations 
are (a) pure a + light, (b) one is a + and the other a~ , and (c) lin _L lin. The overlap region near 
zero delay results in an optical standing wave with a period of 177 nm. The solid curves of (a) and 
(c) are calculations based on independent measurements of the optical pulse chirp and envelope 
with no free fit parameters. The bold curve also incorporates a thermal average of the intial ion 
position at the Doppler cooling limit of n = 40, washing out the fine standing wave fringes. 



orthogonal linear polarizations ("lin _L lin") where neither beam drives Raman transitions 
by itself. 

In order to quantify the population transfer probability, we compared detection his- 
tograms [15] from a fast Raman 7r-pulse to a 7r-pulse applied with microwaves and measured 
a population transfer probability of 99.3%. We also investigated the population transfer 
efficiency of the 7r-pulses by repeating the gate multiple times at the laser repetition rate 
and monitoring the transfer probability while increasing the number of 7r-pulses. Since the 
laser repetition period is very close to a half integer number of hyperfine evolution peri- 
ods (^HF//rep. ~ 104.5), there is some natural spin-echo-type error cancellation for multiple 
pulses for our system. For the a + /a + configuration, we measure a population transfer 
probability contrast of 91% after 25 7r-pulses while for a + jo~ we measure 77%. 

In the overlap region (near zero delay time) in Fig. [31 the two counterpropagating tt/2- 



pulses begin to overlap in time at the position of the atom and form an optical standing 
wave. For the a + /a + and lin _L lin configurations, the two-photon Rabi frequency has 177 
nm period spatial interference fringes. In the overlap regions of Fig. E^a) and (c), Raman 
transitions are sensitive to the motion of the ion since momentum kicks of 2hk are being 
transferred from the optical field. The pulse's spectrum, however, is much wider than the 
trap frequency (100 GHz compared to 500 kHz), so the pulse simultaneously drives many 
motional sidebands (including the carrier transition). For the a + /a~ configuration (Fig. 
Efb)), the overlap region contains linear polarization and does not drive Raman transitions. 

The thin solid curves in Fig. [3] show the results of a numerical solution of the Schrodinger 
equation vs. delay, which oscillate at optical-frequency delays of less than 1 fs in the overlap 
region. In this region, the ion experiences momentum kicks and phase shifts as spin states 
spread out in motional phase space, leaving the interferometer open. Since the final spin 
phase of different motional states depends sensitively on the ion's intial position, thermal 
averaging tends to wash out the fast spatial variation of the transition probability. The thick 
curves show a thermal average over the ion's initial state, assumed to be a thermal state 
at the Doppler cooling limit of n = 40. The thin and thick curves in Fig. [3] are not fit to 
the data as these numerical solutions are fully constrained by the pulse duration and chirp 
measurements shown in Fig. QJb) and [2](a/ . 

Implementation of the fast cooling [lOj and fast entangling gates [111 . Il2| will require 
repeated spin-dependent momentum kicks generated through ultrafast interferometry in the 
strong excitation regime [9(. In this case, the interferometer will be closed and therefore will 
not be sensitive to optical wavelength interference such as the overlap region in Fig. [3]^c). 

We have demonstrated ultrafast single-qubit gates with a mode-locked laser pulse using 
an atomic qubit. For a single trapped ion, the free-evolution of the qubit can be used to 
perform a z rotations, and delaying the pulse arrival time will allow a rotation about an 
arbitrary axis in the x-y plane of the Bloch sphere. As such, the fundamental limit on the 
gate speed is the hyperfine period (analogous to the Larmor precession time), which would 



rese same qubit levels in m Yb + 
so this single-qubit gate can be 



yield a gate time of ~ 100 ps. Previous results with t 

have demonstrated coherence times in excess of 1000 s 

performed in a vanishingly small fraction (< 10~ 13 ) of the coherence time. 
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